L'architecture des équations différentielles linéaires d'ordre supérieur
Une équation différentielle linéaire d'ordre $n$ est caractérisée par sa dérivée d'ordre le plus élevé. Nous définissons la forme générale comme l'Équation (1) :
$$P_0(t) \frac{d^n y}{dt^n} + P_1(t) \frac{d^{n-1} y}{dt^{n-1}} + \dots + P_{n-1}(t) \frac{dy}{dt} + P_n(t)y = G(t)$$ (1)
Pour faciliter l'analyse théorique, nous normalisons souvent cette équation en la divisant par $P_0(t)$, en supposant qu'elle est non nulle sur l'intervalle d'intérêt. Cela donne la Forme standard (Équation 2) :
$$L[y] = \frac{d^n y}{dt^n} + p_1(t) \frac{d^{n-1} y}{dt^{n-1}} + \dots + p_{n-1}(t) \frac{dy}{dt} + p_n(t)y = g(t)$$ (2)
Notation d'opérateur et coefficients constants
La complexité de $n$ dérivées est consolidée en un seul opérateur linéaire $L$. Lorsque les coefficients sont constants ($a_n$), l'expression se simplifie en :
$L[y] = a_0y^{(n)} + a_1y^{(n-1)} + \dots + a_{n-1}y' + a_ny = g(t)$
Cette notation met en évidence que $L$ agit de manière linéaire : $L[c_1y_1 + c_2y_2] = c_1L[y_1] + c_2L[y_2]$. Ce principe garantit que la solution générale est composée d'une solution complémentaire ($y_c$) et une solution particulière ($Y$).
Considérez Figure 4.2.4: Un système à deux ressorts et deux masses avec des masses $m_1, m_2$ et des déplacements $u_1, u_2$. La physique fournit deux équations d'ordre deux couplées. En isolant $u_1$ par substitution, nous obtenons une seule équation d'ordre 4 équation. Pour la résoudre, nous avons besoin de 4 conditions initiales (position et vitesse pour chaque masse) afin de trouver une trajectoire physique unique.
Exemple travaillé : La solution homogène
Trouvez la solution générale de l'équation différentielle : $y''' - y'' - y' + y = 0$
Supposons $y = e^{rt}$. En substituant dans l'EDO, on obtient : $r^3 - r^2 - r + 1 = 0$.
Factorisation par regroupement : $r^2(r - 1) - 1(r - 1) = 0 \implies (r^2 - 1)(r - 1) = 0$.
Cela s'expande en $(r - 1)(r + 1)(r - 1) = (r - 1)^2(r + 1) = 0$.
Les racines sont $r = 1$ (multiplicité 2) et $r = -1$. Comme $r=1$ se répète, nous multiplions le deuxième terme par $t$.
$y_c(t) = c_1e^t + c_2te^t + c_3e^{-t}$